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Accepting the fine tuned cosmological constant hypothesis, we have recently proposed 
that this hypothesis can be tested if the dark matter freeze out occurs at the electroweak 
scale and if one were to measure an anomalous shift in the dark matter relic abundance. 
In this paper, we numerically compute this relic abundance shift in the context of explicit 
singlet extensions of the Standard Model and explore the properties of the phase transition 
which would lead to the observationally most favorable scenario. Through the numerical 
exploration, we explicitly identify a parameter space in a singlet extension of the standard 
model which gives order unity observable effects. We also clarify the notion of a temperature 
dependence in the vacuum energy. 
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1. INTRODUCTION 

The hypothesis that the cosmological constant (CC) energy density today is a result of a tuning 
between UV and IR contributions [U |2] is favored according to some versions of the string landscape 
proposal (see e.g. [3]). Furthermore, this hypothesis has always been the default assumption in 
particle physics model building (see e.g. [HE]). Unfortunately, this conjecture is notoriously difficult 
to test with lab experiments, such as those at colliders. 

One of the predictions of the tuning hypothesis is that there can be an electroweak scale effective 
CC in the early universe if there was a phase transition (PT) at that scale. A well-known reason 
to suspect that there was an electroweak scale PT in the early universe is the thermally supported 
electroweak symmetry restoration phenomenon in the context of the Standard Model (SM) of par- 
ticle physics [HE]. Hence, if lab experiments, such as particle colliders, can eventually measure the 
field content and couplings of the scalar sector at the electroweak scale with sufficient accuracy, then 
one may be able to predict the CC energy density existing at around the time of the PT. Such an 
energy density would interact with gravity to modify the expansion history of the universe. Indeed, 
Kolb and Wolfram [8j were one of the first to state that this computable energy density arising from 
the Standard Model Higgs condensate may have an observationally acceptable yet significant effect 
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in cosmology. 

In a recent paper [9], we proposed that dark matter freeze out can be used to probe PTs, including 
the properties of such a computable CC, through its effect on the expansion rate of the universe 
during freeze out. Such an idea is abstractly very similar to the well known big bang nucleosynthesis 
idea, as well as generic particle probes of cosmology (see e.g. |10H14j ). In particular, if a weakly 
interacting massive particle (WIMP) dark matter candidate is discovered with a mass of the order of 
TeV, then its freeze out dynamics would be sensitive to the value of the CC during the electroweak 
scale PT at a temperature of the order of 100 GeV. Therefore accurate measurements of the dark 
matter and scalar sector properties will, in principle, make possible a lab test of the tuning of the 
CC. More accurately, what is being tested is the absence of self-tuning mechanisms and/or modified 
gravitational theories |15H27| that would eliminate or significantly change the effects of vacuum 
energy during a PT. 

For non-first order PTs, it was found that the shift in the relic abundance due to the CC energy 
density effects is suppressed by Anx/nx = O (g^ 1 ) where is the number of of relativistic degrees 
of freedom contributing to the energy density. For first order PTs, it was found that this fractional 
shift can be generically enhanced by supercooling such that the CC effects can be 0(1) with a 1% 
parameteric tuning. In all cases, the sought after CC "signal" is buried in the dominant "background" 
coming from the adiabatic change in the number of degrees of freedom and possibly the entropy 
release near the time of the dark matter freeze out. The adiabatic change in the number of degrees 
of freedom and the vacuum energy effect tend to increase the relic density today while the entropy 
production effect decreases the relic density. 

The purpose of this paper is to complement the previous short paper [9] in several ways: 

1. Present an explicit singlet extension of the Standard Model (SM) that gives a large super- 
cooling with a first order PT at the electroweak scale. 

2. Clarify the notion of how an effective vacuum energy (which is Lorentz invariant in the flat 
space limit) can depend on temperature (which manifestly breaks Lorentz invariance). 

3. Compare numerical results with analytic results presented in |9]. 

4. Provide relevant technical details that were left out in [9] to aid future research efforts in this 
direction. 

In addition to giving a generic singlet scalar model coupled to a Dirac fermion that gives a significant 
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supercooling, we analyze xSM, i.e. a real singlet coupled to the SM, and identify a parametric region 
in which significant supercooling occurs. As anticipated in [9], an O(10 -2 ) tuning is sufficient to 
induce an 0(1) supercooling effect on the relic abundance. 

The order of presentation is as follows. In Section [2] we review the physics of PTs and focus 
on the myriad ways in which a PT may impact dark matter freeze out. We clarify the notion of 
a temperature dependence of vacuum energy density in this section. In Section [3] we analytically 
compute the fractional shift of the relic abundance 5nx(to) due to an electroweak scale PT in 
the limit in which the PT represents a small perturbation to the usual freeze out. In Section [4] we 
compute the relic abundance deviation in the SM and minimal singlet extensions (both supplemented 
by a generic dark matter which is assumed to play a negligible role in determining the properties of 
the PT). In Section [5] we conclude with a summary and suggestions for future work. An extensive 
set of appendices detail technicalities useful for the material presented in the body of the paper. 

Throughout the paper, we assume a flat Friedmann-Robertson- Walker (FRW) metric, ds 2 = 
dt 2 — a 2 (t)\dx\ 2 , and use the reduced Planck mass M p « 2.4 x 10 18 GeV. 

2. A BRIEF REVIEW OF THE PHYSICS OF PHASE TRANSITIONS 

In this section, we review the physical features that accompany a cosmological PT. Each of these 
features modifies one of the relationships, p ~ T 4 , T ~ a -1 , or (av) = (av) (T), which are assumed 
in the usual freeze out calculation. One of the topics discussed in this section is how to understand 
the thermal dependence of vacuum energy, which a priori is an oxymoron. Readers interested in 
mostly the phenomenology can skip to the next section. 

The standard cosmological model assumes an expanding FRW universe which leads to the tem- 
perature of the relativistic species in the universe decreasing as a function of time except during 
the time periods when entropy is generated. As the temperature decreases, there may exist critical 
temperatures at which the thermodynamic quantities are not analytic as a function of tempera- 
ture and/or the symmetries of the effective Lagrangian governing the dynamical degrees of freedom 
changes. Following the typical convention in the literature, we refer to the passages through these 
critical temperatures as PTs. 

In order to calculate thermodynamic quantities in the system described above, we will use the 
thermal effective potential (see |28] for a review). The thermal effective potential V e s((f) c , T), derived 
from Legendre transforming the partition function coupled to external sources, represents the free 



5 



energy density of the plasma at temperature T dynamically interacting with a homogeneous scalar 
field background 4> c which may affect the masses and interactions of particles in the plasma. A 
local minimum 4> c = v(T) is called the thermal vacuum, and PTs occur near critical temperatures 
T c which will be defined more precisely below. 1 

The critical temperature T c in the case of what is conventionally referred to as a first order PT 
is defined by the existence of two or more degenerate minima <p = v(T c ) existing for the thermal 
effective potential V e d(4>, T c ). In such cases, we refer to the vacuum of the universe just prior to 
the PT as (T) (where the "s" superscript denotes "symmetric" vacuum) whether or not there is 
a symmetry in the thermal effective potential prior to the PT. The vacuum solution after the first 
order PT is referred to as v ( fe )(T) where "b" denotes "broken." A non-first order PT (sometimes 
loosely referred to as a second order PT) is characterized by a single continuous function v (T) before 
and after the PT: i.e. v^ b \T c ) = v^ s \T c ). Even in such situations, it is sometimes useful to define 
v( s \T) to be the vacuum before the PT whenever there is a restored symmetry prior to the PT. The 
quantity v^ s '(T) can then be taken as an order parameter associated with spontaneous symmetry 
breaking. 

The thermal vacua «( S / 6 )(T) can be obtained from summing up thermal tadpole corrections 
obtained from expanding perturbatively about the zero temperature vacua v^ s ^ b \0). Despite the 
suggestive notation of the thermally shifted vacuum v ( s / b \T), the resummation of tadpoles is noth- 
ing more than a reorganization of perturbation theory, and the vacuum energy represented by the 
Lorentz-invariant part of the energy-momentum tensor, is not shifted by the manifestly Lorentz- 
noninvariant thermal tadpoles. Note that symmetry restoration cannot be inferred from the thermal 
tadpole resummation alone since the thermal perturbation theory breaks down when the perturba- 
tions are expanded about the inflection points of the effective potential. 

Let us now establish some more notation for the quantities introduced above. The thermal 
effective potential and ?/ s / 6 )(T) can be used to construct the thermodynamic quantities 

F {s / b \T) = V cS (v (s ^(T),T) (2.1a) 
s (s/b)^ = _d jr{s/b) (2lb) 

pWb)( T ) = jr(s/b) +Ts (s/b) ( 2 lc ) 

representing the free energy density entropy density s, and energy density p in the symmetric and 



1 We will leave out the adjective "thermal" in "thermal vacuum" whenever no confusion should arise. 



broken phases. A typical PT occurs as the universe cools, and the free energy of the broken phase, 
in which the entropy and energy densities are high, drops below the free energy of the symmetric 
phase, in which the entropy and energy densities are low. It will be useful to define the critical 
temperature of the PT T c by 

T^(T C )=T^(T C ), (2.2) 

but note that the PT may not actually occur until a much lower temperature if the symmetric phase 
experiences supercooling. The PT is accompanied by a number of physical features, which we will 
outline in the remainder of this section and which each have an impact on dark matter freeze out. 

The first feature that we would like to discuss is the vacuum energy associated with the PT. 
We assume that the energy density p^ s l b \T) can be partitioned into the energy associated with the 
plasma and the energy associated with the condensate (i.e. the vacuum energy with an effective 
equation of state of —1), and we define the latter as 

p^ b \T) = V eS (v {s/h) ino) (2.3) 

which has an observable consequence when coupled to gravity. This equation is artificial because the 
vacuum energy cannot be rigorously separated from the particle energy with which it is associated 
for most of the states populating the density matrix. Nonetheless, it is useful because it captures 
the CC type of contribution (i.e. negative equation of state contribution) to the energy-momentum 
tensor. 

Note that flat space thermal corrections to the zero temperature effective potential cannot gen- 
erate Lorentz invariant contributions to the energy-momentum tensor because temperature T de- 
pendent quantities are not Lorentz invariant. Since the CC contribution to the energy-momentum 



tensor in the flat space limit is Lorentz invariant, one may wonder whether Eq. (2.3) is valid since 
it implies that thermal tadpoles are contributing to the vacuum energy. Furthermore, the fact that 
the effective vacuum energy takes on a continuum of values while the only non-pert urbatively stable 
vacuum state is at v^(0) (which we will assume to be associated with negligible vacuum energy) 



also leads one to be suspicious of Eq. (2.3). 



To semi-quantitatively resolve this puzzle, one notes that near the time of the PT, there are 
A -H- B processes in equilibrium where A and B schematically correspond to states of the form 
|particles + vacuum energy) and |particles), respectively. These transitions are mediated by non- 
perturbative processes since they are vacuum changing processes. Classically, the plasma (when 
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these transitions are efficient) is approximately described by inhomogeneous solutions in Minkowski 
space. This can easily be characterized by computing for example the thermal two-point function. 

The equation of state for such a plasma in the classical approximation corresponds to neither 
that of quantum expectation values with respect to states A nor B, but is a mixture which from the 
quantum perspective depends on the non-perturbative transition operators as well as the relative 
statistical and/or coherent weighting of A and B type of states. The incoherent aspect of this 



mixture is what the T dependence of Eq. (2.3) reflects. 2 To corroborate this picture, one can 



easily solve classical equations of motion in models with spontaneous symmetry breaking to obtain 
inhomogeneous background field solutions which have an inhomogeneous equation of state. Since 
the Friedmann equation (which is the only gravitational probe we will be concerned in this paper) 
approximately describes the gravitational response to the spatial average of the energy-momentum 
tensor, one can spatially average the energy density and the pressure. This leads to an effectively 
homogeneous energy density and pressure which is approximately the same as that due to particles 
plus a vacuum condensate energy density. The resulting effective vacuum condensate energy density 
is somewhere between V e ff (v^ s \0), 0) and V e g (v^ b \0), 0), justifying the diagnostic quantity defined 



by Eq. (2.3). 



To renormalize the CC, we impose the tuning condition 

pW(T = 0) = 0, (2.4) 

which states that the vacuum energy density today is on the order of the meV 4 dark energy density 
|29[I30| and negligible as compared to the PT scale. Hence, we will refer to Pa/ b \T) as the "effective 
CC energy density." With this normalization, the CC energy density before a PT at scale M will 
typically be 

p£(T>M)~ M \ (2.5) 

which can be measured, in principle, by gravitational probes such as the Hubble expansion rate 
and its impact on dark matter freeze out. 3 Any self-tuning/modified-gravity mechanism which 



2 Note also unlike in flat spacetime, there are IR cutoffs associated with the expansion rate H for a single causal 
domain during the PT and Ho associated with the presently observable universe. The former scale H is also 
associated with one of the scales at wich quasi-equilibrium assumption breaks down. 

3 Although an in depth discussion of UV sensitivity of the CC is beyond the scope of this paper, one should keep 



in mind that using Eq. (2.4 1 as a quantum renormalization condition leads to Eq. (2.5 1 as a prediction only if 
assumptions about analyticity of the effective potential as well as Lorentz invariance structure of the UV cutoff is 
assumed. 
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decouples the vacuum energy or significantly modifies the vacuum energy effect on gravity on a 

(s) 

time scale shorter than that of the expansion scale will have an effective p C c significantly different 



from Eq. (2.5). It would be interesting in future studies to compare various self-tuning/modified 
gravity models which may have non-trivial time dependence in the effective vacuum energy different 
from that in this paper. 

The second PT feature is the decoupling of heavy degrees of freedom which become non- 
relativistic after the PT and cause the number of relativistic species, denoted here as g, to decrease. 
This has two consequences for the dark matter freeze out calculation. First, the energy density of 
the plasma p ~ g T 4 and Hubble expansion rate H ~ ^fp decrease more rapidly than usual after the 
PT. Second, since temperature is related to the FRW scale factor by entropy conservation, which 
gives T ~ <7 _1 / 3 a _1 , the temperature decreases less rapidly than usual after the PT. To estimate the 
magnitude of the effect on dark matter freeze out, consider the SM electroweak PT at T ~ 100 GeV 
and suppose that freeze out occurs at the same temperature. Then during the residual annihilation 
stage of freezeout, which lasts until T ~ 10 GeV, g will decrease by approximately 20% corre- 
sponding to the decoupling of the top, Higgs, and massive gauge bosons. In the usual freeze out 
calculation, changes in g are neglected, because freeze out occurs much later than the electroweak 
PT when g is insensitive to T. When we arrange for the two events to occur at the same scale, g 
decreases significantly and can have a large effect on the relic abundance. 

The third feature is related to the coupling between the PT sector and the rest of the particle 
physics model. As the phase changes at the PT, in general the masses and interactions of particles 
in the plasma can change as well. In particular, it is possible for the scalar field to couple to dark 
matter in such a way that the dark matter's mass and/or annihilation cross section is different in 
the symmetric and broken phases. This scenario, studied by |31| [32], may allow dark matter to 
rethermalize and can have a significant effect on the relic abundance. 

If the PT is of the first order, then it possesses a number of additional features. A first order PT 
can be divided into two stages. The first stage, known as supercooling, occurs while the universe 
remains in the symmetric phase after it has become metastable at T ~ T c . As the temperature 
decreases and the CC energy density remains approximately constant, the total energy density 
can deviate from the standard p ~ T 4 scaling (i.e., first feature above). Supercooling ends when 
it becomes energetically favorable for bubbles of the broken phase to nucleate. Determining the 
temperature Tp T at which bubble nucleation begins requires one to solve for the non-perturbative 
bounce solution and evaluate the decay rate of the metastable phase |33j . During the second stage, 



9 



known as reheating, the expanding bubbles release an energy density 

Ap ex = (Tp T ) - p® (Tp T ) (2.6) 

which is converted into radiation and heats the gas from Tp T before the PT to Tp T > Tp T after the 
PT. We assume that reheating occurs rapidly as compared with the expansion rate of the universe 4 , 
which allows us to treat reheating as an abrupt process at time tpT when a = apT- Reheating is 
accompanied by a non-adiabatic entropy increase. This entropy growth modifies the relationship 
between temperature and the FRW scale factor in such a way that the universe is relatively larger 
for a given temperature. As a result, the dark matter number density undergoes a longer period of 
dilution and the relic abundance can be significantly smaller [35]. Finally, just as massive species can 
adiabatically decouple after the electroweak PT occurs, heavy particles can undergo a non-adiabatic 
decoupling at the time of a first order PT if they abruptly acquire a mass m > Tpx- 



3. AN ANALYTIC ESTIMATE OF THE CHANGE IN THE DARK MATTER 

ABUNDANCE 



In this section, we estimate the change in the dark matter relic abundance due to the presence of 
a PT, and the CC energy density in particular, during freeze out. Our final result is the fractional 



deviation of the relic abundance, denoted 5nx(to) an d given by Eq. (3.44), in which we have 



linearized in the various effects of the PT on freeze out. Although most of the results in this section 
have already been presented in |9], we repeat some of the results for self-containedness as well as 



serving as introduction for more complete results such as Eqs. (3.46) and (3.47). The main point 



of this section is to present a formalism to understand analytically the effects outlined in Sec. |3j 
Throughout the calculation, we will take a as the independent variable and rewrite functions of 



temperature using T = T(a) given by Eq. (D.12). In particular, we will assume that freeze out 



occurs at a temperature Tj = T(df) before the PT at a = apx- Since all of the thermodynamic 
quantities depend on the phase of the system which changes at a = apT, the formulas in this 
section would become unnecessarily obscure if we persisted in writing all the (s/b) superscripts and 
distinguishing the a < q>pt and a > apT cases. Therefore, we introduce the following shorthand. 



4 A third stage, known as phase coexistence, can occur if a large latent heat is released by the expanding bubbles and 
the plasma is reheated to the point where the pressure gradient across the bubble wall vanishes [34] . Subsequently, 
the bubbles expand only insofar as the universe expands, and the PT completes on a time scale t ~ H ■ Typically, 
this stage does not occur during an electroweak-scale PT because the number of relativistic species O (100) is too 
many to allow for sufficient reheating. 
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Whenever a temperature-dependent function F( S / 6 )(T) appears without the (s/b) superscript, the 
intended meaning is 

F( s )(T(a)) a<a PT 



F(a) 



(3.1) 

F( b \T(a)) a>a PT 



In particular, one always has F(a f ) = F^(T f ) since aj < apT by assumption. 

We calculate the thermal relic abundance of dark matter by integrating the thermally averaged 
Boltzmann equation, 

^(A x )=-H(4-4 q2 )- (3.2) 

over the era of residual annihilations from freeze out at a = af until today. Subject to general 
assumptions (see Appendix |B] for more details), we obtain 



nx(to)=[l ' ""^"" M^ (3.3) 



inoo/a/ d\u(a/a f ) , ,al 

jr— {av) a^ 

for the number density of dark matter today at a = a>o and t = to- I n this expression, the quantities 
that will be affected by the PT are the Hubble expansion rate H(a), the thermally averaged cross 
section (av) (a), and the dilution number since the time of the freeze out to today ao/at, which is 
related to T(a). As a fiducial reference value, we also calculate the "usual" relic abundance Uy (to) 
by assuming that the PT does not occur, but instead that the universe remains radiation dominated 
and has the standard scaling relations 

HM ~ a~ 2 , (av) {U) = (av) (T(o)), and ~ a' 1 (3.4) 

throughout freeze out. We define the relic abundance fractional deviation as 

Sn x (t ) = ^§^--1 (3.5) 

and expect this quantity to depend on the way in which H, {av), and ao/ a / deviate from the usual 
freeze out scenario. We will consider each effect in turn. 



Before addressing each of the factors in Eq. (3.3), let us discuss the partitioning of energy. 



The Hubble expansion rate, which appears in Eq. (3.3), is related to the total energy density 
p( s / b )(T). However, we are particularly interested in determining the impact of the effective CC 
on the calculation of dark matter freeze out. Therefore we will assume that the energy can be 
partitioned as 

p 7a (particle degrees of freedom + exotic energy component) . (3-6) 
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In general, the exotic energy component can arise from physics other than the effective CC, such as 
quintessence (e.g. |13 | I14 [ l36Tl40| ) or late-decaying massive particles (e.g. |41H48j ). To maintain a 
minimal degree of generality throughout our analytic estimates (without accumulating distasteful 
notational complication), we will parametrize the exotic energy component as p ex n(a). However, 
since our primary interest is in the case that the exotic energy component represents an effective 
CC, we will write 



/9 CX k(o) = p cc (a) 



(3.7) 



where p^iJ b \T) is defined by Eq. (2.3), and we have used the shorthand Eq. (3.1). The remaining 
energy density can be attributed to relativistic particles in the plasma, which we will denote by 5 



P ( R /b \T) 



P 



(s/b) 



(T) 



Fee 



(T). 



(3.8) 



To connect with a familiar and intuitive notation, we let the functions and g$ be defined implicitly 
by 



s (s/b) (T) 



30 
2tt2 

45 



9 T\t)t" 

g ( s s/b \T)T 3 



(3.9) 
(3.10) 



such that they represent the number of relativistic degrees of freedom at temperature T in either 
the symmetric or broken phase. As shown in Appendix [Cj one must have gs(T) ^ 9e(T) if entropy 
and energy are to be conserved during the time when a species adiabatically decouples. 



Now, we will begin our investigation of the quantities in Eq. (3.3). First, consider the effect on 



the Hubble expansion rate H(a) which is obtained by solving the Friedmann equation. To do so, 
we partition the energy as described above and assume that p ex <C pn(cif) such that we can treat 
the CC energy density as a perturbation. With these assumptions, we obtain 

1 



H(a) 



3M D V 10 



7T^ 



9E{a) 



1 + 



1 



Pcx K(a) 



2^g E (a)T(ar 



(3.11) 
(3.12) 



Contributions from non-relativistic species are Boltzmann suppress ed. D efined in this way, p^ 6 ' includes a term 



proportional to dw^ s//f! ' jdT which arises from the derivative in Eq. (2.16 1. This term represents kinetic energy in 



the scalar field and, strictly speaking, should not be included in pn. Nevertheless, we do not separate out the 
kinetic term, because it is typically negligible. 
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where we have used the shorthand Eq. (3.1). During the PT, we can approximate k(o) as 

A/9 e 



k(o) Q(apT — a) + 0(a — apr) 1 



(3.13) 



where Q(z) is a step function, A/9 CX > is given by Eq. (2.6), and At2(a) is a function which starts 
from K2(apx) = 1 and decreases as fast as 

(^r < 3 "» 

\apT J 

with > 4. If A/j ex = 0, we have a continuous second order transition or a crossover. If Ap cx = p cx , 
then the supercooling is sufficiently strong as to end up with no CC energy just after the PT. The 
step functions represent the fact that the PT occurs with negligible change in the scale factor. 
With this assumption, As and the corresponding change in the temperature become functions of 



Ap ex according to Eq. (D.8) in Appendix [Dj Finally, the Q(apx — a) term in Eq. (3.13) should, 
in general, be multiplied by another smooth function unless there is some symmetry fixing v^ s '(T), 
and consequently p^c(T), to a particular value in the high energy limit. However, we will neglect 
this detail in favor of cleaner notation, since the final result will be approximately unchanged. 

As discussed in Section |2j particle species start becoming non-relativistic after the (electroweak) 
PT which causes gE/s( a ) to decrease. We will parametrize this decrease by focusing on the (non- 
)adiabatic decoupling of (Npt) N fermionic dynamical degrees of freedom and write 



9E/s( a ) = 9E/s(af) ~ K a 
7 



h{a) 



7 



N PT 8(a - a PT ) + -N f(a 



(3.15) 
(3.16) 



where /(a), which rises from to 1, is given by Eq. (EA). Note that in reality, h(a) is a smooth 
complicated function (particularly Nf(a)), but here we are accounting for the change in the number 
of degrees of freedom in a physically suggestive approximation. As we will see below, this effect will 
be one of the dominant "backgrounds" to the "signal" of measuring the effects of the cosmological 
constant. We treat this effect as a perturbation to linear order, and we estimate the Hubble 
expansion rate to be 

T(a) 2 



H{a) 



3M„ 



9E(df) 



1 



1 h(a) 



+ 



1 



Pcx «(a) 



2g E (a f ) 2^g E ( af )T(c 



(3.17) 



Writing T(a) using Eq. (D.12) and linearizing further with respect to small quantities, we have 

l 



H{a)^H (u \c 



1 + 61 - 
2 \a f 



2 11 

K ( a ) + o £ 2©( a ~ a PT) + ^e3i©(a - a PT ) + ^£32/ (a) 

3 6 6 



(3.18) 
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where 



H {U) (a) ee T J— 2X j'^- gE{af) (3.19) 



3M 



P \ a 



and 



ei ee — — = fractional energy of the exotic during freeze out (3.20a) 

h9E{a f )Tf 

( (X \ ^ A .q 

62 = I ) = fractional entropy increase during PT (3.20b) 



i 

INpt 

631 ee -2— — - = fractional decoupling degrees of freedom during PT (3.20c) 

9E{af) 

in 

632 ee — — — r- = fractional decoupling degrees of freedom near freeze out (3.20d) 

9E{a f ) 



where As, denoting the entropy density change at the time of the PT, is given by Eq. (D.8). 
Although H(a) appears to vary discontinuously at a = apx, its continuity is ensured by the con- 
servation of energy. At the PT, the CC energy converts into radiation, which generates an entropy 
but leaves the total energy density fixed (i.e., €2 compensates for the discontinuity of the ei term) 
because the volume remains approximately constant through the duration of the PT. The fact that 
H is boosted by £31 and £32 is intuitive for the following reason. When a particle species becomes 
non-relativistic, the effective equation of state becomes smaller, such that the energy dilutes less, 
which in turn leads to a larger expansion rate for the same scale factor. The term 631 accounts for 
the non-adiabatic change in the number of degrees of freedom during the PT, while the term €32 
accounts for the adiabatic change in the number of degrees of freedom. 

Next, consider the change in the cross section due to the PT. We parametrize this effect as 

(av) = (av) (U) (l - e 4 8(a - opt)) (3.21) 

where 

e 4 EE--^- (3.22) 



and Ao- is the change in (av) due to the PT. Since the derivation of Eq. (3.3) assumes that the 
dark matter is decoupled after Tf, we will assume that £4^0 in order to prevent re-thermalization 



due to an increase in the cross section. Hence, we can evaluate Eq. (3.3) by linearizing in the e's to 
obtain 



nx(to) 



a f\ 3 ( f lnao/af din 0/0/ (av) (U) " ~ X 



a J W H(U)( a ) (a/a f f 



1 + E' 



[a e„ 



(3.23) 
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where 



2 



«(o) - 6(a - a PT )-e 2 



[e 3 i G(a - ap T ) + £32 / (a)] - e 4 G(a - a PT ) 



(3.24) 



implicitly defines the 9 n . Note that the integral is dominated by contributions around In a/ (if = 0. 
On the other hand, the (a//ao) 3 prefactor should be evaluated with all the gs changes accounted 
for, not just the effects around lnx = 0. 

Next, let's consider the effects on the af/ao factor determined by the freeze out condition itself. 
The freeze out temperature Tf can be solved using [8j 



{av)n^(T s ) = C n ^H(T f ) 



Tf 



n£(T) ee g x 



( m x T \ 

v 2vr ) 



3/2 



exp 



T J 



(3.25) 
(3.26) 



where C is an order unity number whose optimum value to reproduce numerical integration is cross 
section dependent (e.g., C 2), gx counts the real dynamical degrees of freedom of the dark matter, 



and mx is the dark matter mass. Evaluating H(Tf) with Eq. (3.18) and assuming freeze out occurs 



before the PT, Eq. (3.25) becomes 

3/2 



(<rv) gx 



exp 



mx \ C mx Tf / tt 2 
~Yf~) ~ 3M P V 10 



9E(a f ) 



1 + 



Although not solvable in closed form, one can linearize in the perturbation again to obtain 



mx_ 
In A 



1 + 



ei 



In A 



+ 0((ln^)" 2 ) 



where 



.4 



g x 3V5M p ^/m x Tf{ 



av) 



exp [20] 



(3.27) 



(3.28) 



(3.29) 



2C^ y /g E {a f ) 

for electroweak mass scales. If we assume that there is only one period of entropy production 
between freeze out and today, and that it occurs at the PT temperature Tpt, we can use entropy 



conservation in the form of Eq. (D.W) to write 



Oo 

a f 



9s(.a f ) \ 1/3 T 
9s(ao) 



_J_ 
To 



1 + 3* 



where Tq is the temperature today. Combining this with Eq. (3.28), we find 



a f 



( a 


(U)\ 


\ a f 





1 ei 1 1 
1+ 2hL4 + 3 62 



(3.30) 



(3.31) 
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where 



ao 
a f 



9s(af)\ 1/3 rn x 1 



gs{ao) J T In A 



(3.32) 



Putting Eq. (3.31) into Eq. (3.23) results in 
-3 



nx(t Q ) 



ao 
a f 



n 3ei 1 
^^lnA" 62 



Ma /a f \W) d l na / af ( av )(V) 

H(U)(a) {ajaff 



where the endpoint contribution to the integral has been written as 



E, = 



a J i_ ^2 

2 In A 3 



+ f M 



(£0 



(3.33) 
(3.34) 



i^)(a ) (oo/o/l^)) 3 ' 

The term is negligible because of the volume dilution factor in its denominator. Linearizing the 
small factors gives 

3ei 1 



n x {to) ~ rJp(t ) 



1 



2 In ,4 



£2 - F u 1 (9 n e n 



(3.35) 



where 



n x (*o) = — 



P^j (a/a/) 3 



F u 1 = usual computation of relic abundance 



a 



(3.36) 



(3.37) 



(3.38) 



, ffd0(a) (a/a/) 3 

In particular, if we assume an s-wave cross section (i.e., constant (cry)), we can express 6\ explicitly 

as 



F- l ~e 



;i + 3.y) 



i 



(3.39) 



n rf - 3 \ /°ex 

where we have expanded in 5 = a pt /a / — 1 > which represents the delay between freeze out and 
the PT. The first term in square brackets comes from integrating the CC energy density from a/ to 
apT, and the second term comes from integrating the decreasing CC energy density after the PT. 
This equation shows that if ~ \ (large supercooling) there is a suppression of the ei effect by 
a factor of order 5. Although we have linearized in 5 along with ej, terms of the form 5 are not 
higher order. The expansion in ei reflects the fact that we treat the PT as a perturbation, whereas 
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the expansion in 5 is performed merely to simplify the expressions. With the same assumptions, we 
can evaluate the other F~ l 9 n terms: 

2 



(1-5) 
1-6) 



1 fH°o/*f\ {U) ) dlna/a f 
(a/a/) 2 



6 JO 
(1-5). 



/(«) 



(3.40) 
(3.41) 

(3.42) 
(3.43) 



Hence, for s-wave cross sections, the change in the relic abundance due to small changes made by 
the PT can be expressed as 



5n x (t ) = ci ei + c 2 e 2 + c 3i e 3 i + c 32 e 32 + c 4 e 4 



(3.44) 



where 



2 \ n-3 V Pex 
C2 = -1(1 + 25) 

c 3 i = i(l-5) 

1 /■M«o/"/l (C7) ) din a/a f , 

C32 = « / / / N 2 /(°) 

6 Jo (a/a f ) z 
C4 = 1 — 5 . 



3 1 



2 In A 



(3.45a) 

(3.45b) 
(3.45c) 

(3.45d) 

(3.45e) 



The key point of Eq. (3.44) is that despite the "background" represented by e n -^i, the "signal" 
contained in e\ can be "measured" and represents a prediction of the hypothesis of a tuned CC. It is 
a tuned but striking statement, nonetheless. Since this term is central to the rest of our calculation, 
we have reproduced the so called "CC effect" term here as 

'(l + 5f-l (l + 5f 



1 



t9E(a f )n T ) (! + *)< \2 



+ 



nd ~ 3 



1 



3 

2 In 



(3.46) 



without linearizing in 5. We also write the so called "entropy effect" as 



5+\ 



As 



c 2 e 2 



(3.47) 



6 + l2 -£gs(ap T )n T 

Note finally that we can obtain a smooth non-first order PT by taking the limit Ap ex = e 2 = e 3 i 
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One should remember that all of the analysis has assumed that the entropy released from the 
PT (in the case of a first order PT) did not reheat the system to the point that the dark matter 
rethermalized after freeze out, i.e., Tp T < Tf. This provides a lower bound on 5 for a given Ap cx , 
which can be expressed as 

1 e 3 i + 1 2 — < $ (3.48) 

4 31 4|J^( a/) (T- T )4~ 



by using Eq. (D.3) and assuming that f(apr) is negligible. 

Note also that the range and independence of {ei,S} that is achievable by choosing a beyond 
the SM Lagrangian is not easy to compute nor to generalize. For example, suppose we want to 
increase 5 while keeping e\ fixed. To increase 5, we increase apT more than a/. Since a/ is mostly 
determined by the mass of the dark matter mx while app is determined in part by the competition 
between the thermal mass support and scalar field mass at the field origin, we can keep a/ fixed 
and increase app by decreasing the scalar field mass competing with the thermal support. This, 
however, typically changes the fractional entropy increase £2 during the PT. Furthermore, this will 
change the index rid (defined in Eq. ( 3.14[ ) ) which depends partly on the flatness of the non-thermal 
part of the scalar potential. Indeed, we see that if this rid can be engineered to be as close to 3 
as possible (i.e. a flat potential with no thermal particles decoupling), then the e\ signal can be 
enhanced. One also sees that in the case of a first order PT, the prediction for the effect of the 
cosmological constant (i.e., the ei piece) depends on Ap ex and 5, both of which depend on knowing 
exactly when the PT occurs. As described in Sec. [2] an accurate computation of this will require a 
non-perturbative numerical treatment. Hence, the first order PT situation, which can give a larger 
CC dependent signal, presents an interesting computational challenge of its own. 



4. ILLUSTRATIVE MODELS 



In this section, we present numerical calculations of Snx(to) for various models. This section 
represents one of the key features of the paper that distinguish it from [9], as discussed in the 
introduction. For each model we specify the parameters of the scalar sector, which appear in the 
thermal effective potential V e s((f) c ,T), and the parameters of the dark matter sector, mj, gx- and 
(crv). We then calculate the relic abundance shift using the methods of Section [3j Most of the 
numerical results have not been reported previously, and the model dependent analysis of a real 
singlet extension of the standard model is entirely new. 
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Figure 1: The energy density at the SM PT, p(T) = p cc (T) + pn(T), relative to the energy density when 
the entire SM is relativistic, p r (106.75) = §^106.75 T 4 . Just before the PT at T > 150 GeV, the energy 
density grows relative to p r due to the temperature independent CC contribution, p cc (T > Tpy) rs const. 
After the PT, the top, bottom, Higgs, and massive gauge boson adiabatically decouple causing pj p r to drop 
below one. This adiabatic decoupling is the dominant feature of the SM PT that is relevant for freeze out. 



4.1. Standard Model with Dark Matter 

We calculate here the relic abundance deviation due to the SM electroweak PT. The qualitative 
results were already given in [9j. The numerical details that we discuss in this section can be 
summarized as 5nx(to) = O (10~ 3 — 10~ 2 ) with the CC contributing c\e\ = O (l0~ 4 — 10~ 3 ) . 
With mh = 115 GeV, the largest CC effect occurs for mx ~ 4.2 TeV where c\e\ ~ 9.5 x 10 -4 . 
Our results are summarized in Figures [2] and [3] In this section, we first discuss these figures and 
then extend the analytic estimate of Section [3] now in the context of a concrete model, to obtain 



Eq. (4.7), which lets us motivate extensions of the SM that achieve larger 5nx(to)- Some of the 
qualitative discussion of [9j is reproduced for completeness. 

In Appendix [F] we compute the SM thermal effective potential V e g(h c ,T) through one-loop 



order 6 , where h(x) = </2\H'H\ is the radial component of the Higgs field and h c = (h(x)). It 



is important to point out that the renormalization conditions, given by Eq. (F.5), are chosen such 



6 It is well known that the one-loop approximation breaks down at the temperature of the SM electroweak PT [49 , 
and that accurate results require lattice calculations 50-52 . However, since the CC contribution already represents 
perturbative correction to dark matter freeze out, we will neglect higher-order corrections to the PT physics and 
simply apply the mean field approximation described in Section [2] 
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Figure 2: Fractional deviation of the relic abundance due to the SM electroweak PT. The numerical 



calculation is represented by the black curve, the analytic estimate Eq. (3.44 1 by the red curve, the CC effect 



(ci€i term) by the green curve, and the adiabatic decoupling effect (032632 term) by the blue curve. The 
right axis shows the values of the c\t\ curve only, and the left axis shows the values of the three other curves. 



that V e g(h c , 0) has a minimum at v = 246 GeV where the curvature is and, most importantly 
the CC is tuned by requiring V e g(v,0) = 0. 

Before discussing the numerical results, it is useful to recall from Section [3] that for a non- 
first order PT, freeze out is only affected by modifications to the relations H(T) oc \J p(T) oc T 2 
and T 3 oc 1 a 3 ~ a 3 . These modifications arise when the energy partitioning deviates from 
radiation domination and the number of relativistic degrees of freedom deviates from a constant 
value. These deviations can be visualized in Figure]!] where we plot p(T) normalized by /9 r (106.75) = 
(vt 2 /30)(106.75) T 4 , the energy density of the SM as if all particles were relativistic. We have taken 
nih = 115 GeV which gives a PT at Tpx ~ 148 GeV. As the temperature decreases toward Tpt 
from above, p/p r grows to approximately 1.006 due to the presence of the additional CC energy 
density (i.e., Au 4 /(106.75 T PT ) m 10~ 3 ). Below T PT the massive species decouple, the plasma loses 
about twenty relativistic degrees of freedom, and p/p r decreases to approximately 0.8. This figure 
illustrates that the adiabatic decoupling has an effect on p which is two orders of magnitude larger 
than that from the CC. Therefore, we expect that the Standard Model electroweak effective CC will 
have a subdominant effect on the relic abundance as well. 

The fractional shift Snx(to) is calculated using the perturbative, analytic expressions in Section 
plas well as by solving the Boltzmann equation numerically. In Figure^] we have plotted 8nx(to) by 
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Figure 3: Fractional deviation of the relic abundance due to the SM electroweak phase transition, plotted 
against the WIMP mass and S. The colors and axes are the same as in Figure [2] 



varying mh and fixing mx = 6 TeV, gx = 2, and (av) = 2.33 x 10 - 39 cm 2 . As seen in the figure, the 
PT causes an O (l0~ 3 — 10 -2 ) fractional increase in the relic abundance. We have chosen the DM 
mass to be 6 TeV such that freeze out and the PT coincide at T ~ 303 GeV for = 300 GeV. 
For smaller mh, the PT is delayed with respect to freezeout. The analytic estimate, given by 



Eq. (3.44), only receives contributions from the CC effect {c\e\ term) and the adiabatic decoupling 
effect (032632 term), because the PT is not first order. As we anticipated in the discussion of the 
preceeding paragraph, the £32 term dominates. The analytic formula consistently underestimates 
the numerical calculation by 2 — 3%, and moreover, in the large limit where 5 ~ 0, the deviation 
grows to approximately 4.5%. Both of these features can be traced back to approximations we 



have made in the analytic estimate. The first is associated with the approximation Eq. (B.5), 
which assumes the number density per comoving volume decreases significantly due to residual 
annihilations and introduces an O (Tf /mx) ^ 5% error at all m^. The second is associated with 
neglecting the equilibrium term in Eq. (3.2), which is not negligible at the start of the residual 
annihilation era. The scaling with m^ also has a simple, intuitive explanation. One can understand 
why 5nx(to) is small at small mh, because in this limit the PT occurs too late and becomes 
decoupled from freeze out. Considering the opposite limit, one may wonder if 5nx(to) continues 
to increase for m^ > 300 GeV where 5 < 0. For 5 < the PT occurs before freeze out, as in the 
usual cosmology, and one would naively expect 8nx(to) = 0. Nevertheless, 5nx(to) does continue 
to grow because of the way we have defined Tl x ■ To calculate the usual relic abundance we 
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assume that there are 106.75 relativistic species at freeze out. If the PT occurs much earlier, the 
number of relativistic species at freeze out will be significantly less than 106.75 and 5nx(to) will be 
non-zero. The CC contribution grows monotonically with decreasing m^, since in this limit the PT 
temperature decreases and c\€\ ~ p ex /Tp T . 

In Figure [3] we plot the relic abundance shift by fixing nih = 115 GeV and varying mx- At 
large mx, freeze out occurs well before the PT, the two events decouple, and the relic abundance 
shift is small. At small mx < 2.8 TeV, freeze out occurs after the PT, and the analytic estimate 
fails. The CC effect c\e\ has a maximum of approximately 10 -3 at b max ~ 0.5. For 5 > 5 max the 
factor ex, given by Eq. (3.20a), is small because Tf in the denominator is large. For 5 < 5 max the 



factor ci, given by Eq. (3.45a), is small because the CC is only present over a short time during 
WIMP residual annihilations. The presence of this maximum suggests that c\t\ will typically be 
more sensitive to variations in the parameters of the scalar sector (e.g., m^) than in variations of 
the DM sector (e.g., mx)- With this in mind, we will focus the remainder of our discussion on 
determining the conditions that a scalar potential must satisfy to maximize c\e\. 

We will now extend the estimates of Section [3] in order to understand Figure [2] through a simple 



analytic approximation. We focus on the CC contribution to 5nx(to), given by Eq. ( |3.46 ), which is 

$n x (t ) 3 ciei ~ P ™ (4.1) 
10 g E Tp T 

up to multiplication by an O (1) function of 5. The factor of g E ~ 106.75 represents the SM 
relativistic degrees of freedom before the PT. If we assume that before the PT, the SM particles 
are light with respect to the temperature, then we can approximate V e s using the so-called high- 
temperature approximation 



V eS (hc, T) ~y~ (h 2 c -v 2 f + cT 2 h 2 c . (4.2) 



Here we have defined 



\ eff = -j (V r eff (0,0)-V r eff(v,0)) (4.3) 

to be the one-loop effective self-coupling and 2cT 2 is the thermal mass acquired by Higgs par- 
ticles passing through the plasma. In the SM and subject to our renormalization scheme, these 
dimensionless numbers are \ e j t ~ \sm an d c ~ C SM where 

csm = (^ m t + 6m fe + 6 m w + 3m z + \ m2 h S j ~ 0,18 (4.4a) 

Asm = ^ + 12g ^4 (48m* + 48m£ _ 2 4m 4 w - \2m\ - (15 + log4)m£) « 0.12 (4.4b) 
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for rrih ~ 115 GeV. The PT occurs at a temperature Tpt where <9^ c Vrff(0, Tpt) = 0. Solving for 
this temperature one obtains 

cT 2 PT = ^fv 2 . (4.5) 

Before the PT, the CC energy density is 

p ex = V eS (0,0) = ^v 4 (4.6) 



and we can estimate the deviation in the relic abundance using Eq. (4.1) to be 

1 1 c 2 ,A ^ 

ciei~— — v • (4.7) 

10 9e A e // 

For natural couplings one expects c 2 /A e // ~ O (1) (e.g., Cg M /XsM ~ 0.28) and finds c\e\ ~ 
1/(10^) ~ 10~ 3 . Recalling also that A e ff ~ m|, one sees that this estimate agrees well with both 
the magnitude and scaling shown in Figure [2] Note that in the X e fj —> limit, we find that both p ex 
and Tpt approach zero, but the ratio p ex /Tp T becomes large. This simple approximation suggests 
that the region of parameter space that maximizes the CC contribution to 5nx(to) will have low 
temperature PTs. This is evident in Figure [2] because the CC effect grows at low where the PT 
temperature is low. Hence we will next consider a model in which a scalar singlet coupled to the 
Higgs is introduced to lower the PT temperature. 



4.2. SM Singet Extension with Z 2 

In this section, we briefly discuss an extension of the Standard Model in which the presence of 
an additional scalar field modifies the electroweak PT dynamics. However, we show that the dark 
matter relic abundance is not significantly enhanced, and we argue that we should consider models 
with first order PTs. Consider an extension of the SM in which a real, singlet, scalar field s(x) 
is coupled to the Higgs h(x) through interactions which respect the Z2 symmetry s — > —s. The 
renormalized potential for this theory can be written as 

U({h, s}) =g {I? - , 2 ) 2 + + + f * 2 (h* - , 2 ) (4.8) 

such that d h U({v,0}) = 0, d%U({v,0}) = m\ and d^U({v,0}) = m 2 s . We require 

2 

m 2 s - a 2 v 2 > and 2a 2 + 64 + 7A > (4.9) 



23 



to ensure (s) = 0. This model, known as the Z2XSM, has been previously studied in order to 
determine the viability of s as a dark matter candidate |53H58| . We will not restrict ourselves to 
this scenario, but instead treat the dark matter as a separate sector. The role of s is simply to 
modify the PT dynamics. 

Since this model possesses a greater parametric freedom than the SM, we can attempt to verify 



the relationship Eq. (4.7), derived in the previous section, which relates c\€\ ~ c /A e //. This is 
accomplished by first mapping the parameters of the Z2XSM to c and \ e ffi and second by performing 
a parameter scan while calculating ciei. We obtain c and X e ff by calculating the thermal effective 
potential as described in the previous section (see also Appendix [F]) . If we assume that the quanta 
of s(x) are light with respect to the temperature, we can then extract c and X e ff by matching the 



effective potential to Eq. (4.2). Doing so yields the expressions 



c = c 5 M + g (4.10) 

2 / 2 

a 2 I a 2 V 



m S 



Ae// = A5M-^V(^) (4.11) 



^(ar) = 3---2(l-- + ^)log[l-z] (4.12) 



X \ X X 2 



where the terms containing 02 arise from 1-loop diagrams with an s-particle in the loop, and the 



function varies from ^(0) = to ^(1) = 1. As a result of the minus sign in Eq. (4.11 ), there is an 
upper bound 0,2 < 5 given by the constraint A e // > 0. Now we can see the impact of the singlet field 
on the PT. For 02 > 0, the parameter c is slightly larger and \ e ff is slightly smaller than in the SM. 



Recall that the PT temperature, given by Eq. (4.5), scales like Tp T ~ \ e ff/c. Hence, the singlet 



field lowers the PT temperature and makes the CC energy density relatively more significant, which 
causes the relic abundance shift to be greater. 

To verify these analytic arguments, we calculate the PT temperature and c\e\ numerically over a 
region of the theory space. We allow m? and 02 to vary in the ranges m\ G [(50 GeV) 2 , (300 GeV) 2 ] 
and a.2 G [—0.1,4.0], and we fix 64 = 0.25 and m 2 = (500 GeV) 2 . The range for is chosen to 
prevent the Higgs from becoming unacceptably light 7 , while the range for 02 is chosen to satisfy 



Eq. (4.9) and to avoid the unitarity bound. We map and ai to c and \ e ff using Eqs. (4.10) 



and ( 4.11| ). In Figure |4j we have plotted the contribution to 5nx(to) from the CC effect (ciei) 



over the c 2 / X e f f—rnh plane. This figure shows that the CC effect grows with increasing c 2 /A e // and 



7 Mixing with the singlet does not significantly reduce the LEP Higgs search bound |54| . Moreover, for small mh the 
electroweak breaking minimum may become metastable [531 159] . and the PT becomes first order |60| . Nevertheless, 
we have allowed mh to be as small as 50 GeV to illustrate the parametric dependence of the CC effect. 
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Figure 4: The CC contribution to Snx(to), given by the c\t\ term of Eq. ( 3.44 1 , plotted over the c 2 / X e j f-rrih 
plane for three values of mx ■ The black line represents the SM (da = 0) . 



decreasing rrih, as we anticipated in Eq. (4.7). The largest value of c\E\ is approximately 1.3 x 10 3 , 
which is only about 40% larger than in the SM. The insignificant enhancement can be understood 



by observing that although a% > tends to decrease c, given by Eq. (4.10), its contribution is 
suppressed by a factor of 24. Since csm ~ 0.18 we run into the unitarity bound on ai before it 
contributes significantly to c. If we were to add N light singlet fields instead of one, the contribution 
to c would be Na>2/24, which can be order one even for small 02- We have not take this approach 
here because the ./V additional relativistic degrees of freedom would have a larger effect on the relic 
abundance by increasing the energy density of radiation than through the CC. We have also plotted 
c%e\ for three different values of the WIMP mass from 4 to 8 TeV. This narrow range of viable 
parameters illustrates the tuning that is required to ensure that the PT and freeze out occur at 
the same time. If the WIMP mass is too large, freeze out occurs too long before the PT when the 
CC energy density was subdominant to the energy density of the plasma. As the WIMP mass is 
lowered, the delay between freeze out and the PT decreases and c\€\ grows. If the WIMP mass 
is too small, freeze out occurs after the PT when the CC energy density has been converted into 
radiation. This is the case in the > 200 region of the mx = 4 TeV plot. 

The examples of the SM and the Z2XSM demonstrate that it is challenging to obtain c\e\ larger 



than O (10 3 ). Our discussion at the end of Section 4.1 and simple dimensional analysis illustrate 



why this is the case. In that calculation we obtained Eq. (4.7) which can be written schematically 
as c\t\ ~ Pex/Tp T ~ c 2 /X e ff. Note that the mass scale v, which controls both p ex and Tpx, cancels 



out in the ratio p ex /Tp T . In light of Eq. (4.7) we propose that the CC effect can be enhanced by 
working in a model that has multiple mass scales if there exists a hierarchy between them. We will 
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explore different applications in the remainder of this section. 



4.3. Generic Single Scalar Model 

In this section we calculate the CC contribution to the relic abundance shift in a generic single 
scalar model. Although extensions of the Standard Model typically contain multiple scalar degrees 
of freedom related by symmetries, the thermal dynamics (supercooling and reheating) of a symmetry 
breaking PT can often be modeled by a single scalar degree of freedom which does not display the 
symmetries of the full theory [9l [61]. With this motivation in mind, we consider the theory of a 
real scalar field <p(x) coupled to N Dirac fields ipi(x). The scalar field will experience a first order 
PT during which dark matter freezes out, and the light fermions will compose the hot thermal 
bath. Using this construction, we will be able to calculate the CC effect, which is related to the 
non-thermal energy density and the amount of supercooling, but we cannot estimate the entropy 
and decoupling effects since these depends on how ip is coupled to the full theory. Therefore, in this 
section we assume no decoupling occurs near the time of the PT and that the number of relativistic 
species is fixed to Qe/s ~ 106.75, the relativistic SM background. Let the action be given by 

%>] = f d A x | i (dp) 2 -U{tp)-Y^i [i$ -mi- hup) ipi + £ ct | (4. 13) 

where 

%)= fe + ^V-V + ^ 4 (4.14) 

is the renormalized potential and £ ct is the counterterm Lagrangian. Note that we have eliminated 
the tadpole term in U(ip) by defining the origin in field space appropriately, but there is still a 
counterterm for the tadpole in C c t- As discussed in Section |4.2| we expect that there will be a 
greater impact on the dark matter relic abundance if freeze out occurs during a first order PT with 
large supercooling. Hence, we would like to understand what region of parameter space yields a PT 
of this kind. In particular, we expect that large supercooling can be obtained if the theory S[p] 
possesses two vacua, which will correspond to the low- and high-temperature phases, and that the 
vacua are separated by a barrier. 

We can determine the vacuum structure by identifying the minima of the effective potential, 
which is calculated in Appendix [F] Provided that the non-thermal radiative corrections are negli- 
gible, the effective potential can be approximated as V e s(p c ,T = 0) ~ U(ip c ). It is convenient to 
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Figure 5: An illustration of the ao dependance of the potential given by Eq. ( |4. 14 1. The curves represent 
ao = —2 (red), ao = (green), ao — 0.5 (blue), ao — 1 (purple), and ao = 9/8 (black). 



eliminate M 2 for the dimensionless quantity ao = AM 2 /2£ 2 while assuming \£ ^ 0. We now see 
that the parameter ao controls the shape of the potential U(ip): for ao = 1> the potential has two 



degenerate minima at ip c = and ip c 



=1 where 



3£ 
2A 



1 + A/1 



!■) 



a 



(4.15) 



for ao > 1, (f c = is the global minimum; for < ao < 1, tp c = v is the global minimum; and 
for ao < 0, (f c = becomes a maximum (see also Figure [5]). Therefore, provided that we take 
< ao < 1, the theory possesses a metastable vacuum in which ip c ~ and a stable vacuum in 
which ip c ~ v. In the stable vacuum, we impose the tuning condition V e s(v, 0) = to solve for 

3/2" 



Pc 



8A 3 



27- 36 a + 8a5 + 27 1 



+ o(n) 



(4.16) 



which represents the CC energy density prior to the PT. Finally, the barrier separating the two 
vacua has a "height" 



[/(barrier) - U{Q) 



1 + O (a ) 



(4.17) 



relative to the metastable vacuum. Due to the factor of a^, the barrier vanishes rapidly as ao 
approaches zero. This is illustrated by the ao = 0.5 curve of Figure [5] in which the barrier is already 
almost imperceptible to the eye. 
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Having established that this theory admits two vacua, we will study the PT using the thermal 
effective potential. Although the numerical calculations use the full effective potential, we can gain 
some intuition by making the high temperature approximation. We assume that the ^-particles 
are light with respect to the temperature of the thermal bath, mf <C T 2 , and that the (^-particles 
are heavy. In this limit, then the one-loop thermal effective potential may be approximated by the 
high temperature expansion 

V cS (p c , T) « U(<p e ) + cTVc + O {m 2 /T 2 ) + O (H) (4.18) 

where c ~ YliLi /12 is related to the couplings between ip and ipi. Just as we introduced ao to 
reparametrize V e s((p c ,0), we can now introduce 

Ac 



a(T) = a l^l + JT^T 2 ) > "o (4.19) 

to parameterize V c g(ip c ,T). This definition is particularly convenient, because now Figure [5] also 
illustrates the temperature dependence of V e g (</? C ; T) (up to (^-independent terms) if one replaces «o 
with ot(T). We obtain the expectation values of ip in the "symmetric" and "broken" phases, v^ s \T) 
and i)^(T), by solving (d / d(p c )V e ff((p c ,T) = subject to the boundary conditions v^ b \0) = v 
and t/ s )(0) = 0. We use the terms "symmetric" and "broken," eventhough S[<p] does not display a 
symmetry in order to connect with the notation of Section [2] 

Provided that this model experiences a first order PT, the CC's effect on the relic abundance 
will depend sensitively on the amount of supercooling at the PT. This is seen by the factor of 



(J/) 4 « (Tp T ) in Eq. ( 3.20a ). Therefore, we will begin by investigating the parametric dependence 
of the amount of supercooling, and we will see that it has an interesting dependance on the parameter 
oiq. The supercooling stage begins when the temperature drops below 



'J>*S\I?-^, (4.20) 



defined by Eq. (2.2), or equivalently when a(T c ) = 1. During supercooling, the universe remains 
in the metastable, symmetric phase until bubbles of the broken phase begin to nucleate. Bubble 
nucleation is a non-perturbative process [62], and it occurs at a rate per unit volume which carries 
the standard exponential suppression V ~ T 4 exp [-S^ /T] , where S( 3 ) (T) is the action of the 0(3) 



symmetric bounce [63-65J. Provided that V e a(ip c ,T) can be expressed in the form of Eq. (4.18) 
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then is well approximated by the empirical formula [66 



^ 13,|(=f/(o. (4.21, 



T T \X 

a ( 2.4 0.26 \ 

with a = a(T). Bubbles form rapidly once the bubble nucleation rate averaged over a Hubble 
volume r// -3 is comparable to the Hubble expansion rate H ~ T 2 /M p . For an electroweak scale 
PT, this equality occurs when S'W/T drops below approximately 140 |67[ 168] , Therefore, we can 
determine the amount of supercooling by solving /T ~ 140 for T = Tp T and comparing this 
temperature with T c . 

Considerations of the equation /T ~ 140 demonstrate that the nature of the PT is strongly 
dependent upon the vacuum structure of the theory as parametrized by ao- We will discuss the two 
cases ao > and ao < separately. For ao > 0, the vacuum with ip c = remains metastable as 
T — > 0. This implies that Tp T can be arbitrarily low, and in this limit of large supercooling the CC 
effect may be arbitrarily large. However, in this case the barrier in F e |j(^ Cl T) persists as T — > 0, 
and it is possible that the PT does not occur at any temperature, but instead that the universe 
becomes trapped in the metastable vacuum. This follows from the observation that for ao > 0, 
/T has a minimum at finite T: at low temperatures S^> jT grows due to the explicit factor of 
T in the denominator, and at high temperatures /(a) diverges as T approaches T c and a — > 1. For 
ao ^ 1 the inequality /T < 140 is not satisfied at any temperature, and the PT does not occur 8 . 
Therefore, if we require that the PT must occur via thermal bubble nucleation, we obtain an upper 
bound on ao- On the other hand, for the case ao < 0, the PT necessarily occurs at a temperature 
Tp T > 0, since the symmetric phase becomes perturbatively unstable at low temperatures. This 
latter case has the drawback that supercooling cannot last an arbitrarily long time. 

Assuming that the PT does occur, we can quantify the amount of supercooling using 

Ssc = 1 - ^ , (4.23) 
which takes values between and 1. Parametrizing the temperature dependance with 5sCi we can 

8 At least, the PT does not occur as a thermal process, although it may still occur as a quantum tunneling process 
|64| . However, since quantum tunneling typically proceeds on a longer time scale, the universe could enter an 
inflationary phase, which leads to a cosmological history that deviates significantly from the perturbations we 
consider in Section [3] 
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A/V7 = 0.20 
A/ Vc r = 1-00 
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1.0 



Figure 6: The amount by which the PT temperature drops below the critical temperature, quantified by 
Ssc , is plotted against the parameter ao which controls the height of the barrier. The curves represent 
X/yfc = 0.04 (purple), 0.20 (green), 1.00 (blue), and 5.00 (red). The square indicates the especially tuned 



parameter set given by Eq. (4.281 



rewrite Eq. (4.21) as 



5(3) 



13.7 



A 



a 



3/2 



/(«) 



(4.24) 



y/cj y/1 - a 1 - 5sc 
a + (l-a )(l-5sc) 2 , (4.25) 
which is now only a function of ao, X/\/c, and 5sc- Of course, this expression is approximate, since 



PT 

a 



we assumed V e g took the form of Eq. (4.18), but it suggests that the amount of supercooling will 
depend most sensitively on ao and A / y/c. Now using the full thermal effective potential, we impose 
/T\ T - = 140 and solve for 5sc, which we have plotted in Figure |6j for various parameter sets: 

£ = 5 GeV A = {0.004, 0.02, 0.10, 0.50} 
N = 1 m = 10 GeV h = 0.346 c w 0.01 . (4.26) 

The supercooling grows with increasing ao and decreasing \/\/c as the barrier height and bounce 
action are made larger. The amount of supercooling is typically 5sc = O (0.5) which implies Tp T = 
O (T c /2). Above a finite value of ao (indicated by a dot) the barrier becomes insurmountably large, 
and the universe becomes trapped in the metastable vaccum. The largest amount of supercooling is 



achieved for X/\/c <C 1 and < ao "C !■ In this parameter regime the CC is large (see Eq. (4.16)), 



and the metastable vacuum is separated from the true vacuum by small barrier (see Eq. (4.17)). 
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Figure 7: The CC effect on the relic abundance c\t\ plotted against cto for mx = 17 TeV (solid), 6 TeV 
(dashed), and 0.3 TeV (dotted). For the contours which are absent, freeze out occurs after the PT when the 
CC is not EW-scale. 



Having come to understand the parametric dependance of the amount of supercooling as <p 
experiences a first order PT, we turn our attention back to calculating the impact of such a PT on 
dark matter freeze out. Using Eq. ( 3.46[ ) we calculate the effect of the CC on the relic abundance 



shift and present the results in Figure [7] We have chosen the same parameters as indicated in 



Eq. (4.26) and have taken 

m x = {0.3, 6.0, 17} TeV g x = 2 (av) = 2.33 x lO^cm" 2 (4.27) 

as well. The figure illustrates that is possible to achieve c\e\ = O (0.01) in the tuned parametric 
regime where A/ 'y/c is small and ao approaches its maximally allowed value. Some of the curves 
are absent for the smaller WIMP masses. This occurs because as mx is lowered, the temperature 
of freeze out decreases as well. In the case that A/ 'y/c is small and the PT temperature is high (see 



Eq. (4.20)), freeze out will occur after the PT for small mx- This statement about the relative 
times of freeze out and the PT also explains why c\e\ is insensitive to ao for certain parameter sets 
(e.g., X/y/c = 1, mx = 17 TeV) and very sensitive for others (e.g., X/y/c = 5, mx = 0.3 TeV). In 
the first case, freeze out occurs long before the PT while in the latter case, freeze out occurs just 
before and during the PT and there is a large impact on the relic abundance. 

To conclude this section, we present a particular tuned parameter set which yields c\e\ = O (1). 
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Suppose that we have only one fermion ip and the parameters of S[(p] are given by 

A = 5.4xl0" 4 h = 0.1 

£ = 0.27 GeV M 2 = (1.89 GeV) 2 m = 10 GeV (4.28) 

which leads to 

v ps 1497 GeV a ps 0.007 

c ss 8.3 x 10~ 4 -4= « 0.018 (4.29) 

and PT temperatures 

T c « 374 GeV T pT ps 16 GeV <5 5C ~ 0.96 . (4.30) 

This parameter set is represented on Figure [6] by a square marker. In the dark matter sector we 
take 

m x = 600 GeV g x = 2 (<rw) = 2.33 x 10~ 39 cm~ 2 (4.31) 

such that 

Tf ?s 34 GeV and 5 ps 1.12 . (4.32) 

Using these values we can estimate the CC effect as 

ciei ps 6.1 (4.33) 

Note that the potential obtained with these parameters has a very shallow metastable vacuum at 
cp c ps 0, separated from the global vacuum at ip c ps u by a very small barrier. 

4.4. Singlet Extension with First Order PT 



In this section, we consider a generalization of the SM extension studied in Section 4.2 in which 
we do not impose a 7L<i symmetry on the singlet field s(x). This leads to model known as the xSM 
|53[ I54j . The xSM admits a first order electroweak PT |61[ [69], and we seek to compute the effect 



on the relic abundance due to the effective CC at the PT. As discussed in Section 4.3 the CC 
effect grows with the duration of supercooling. With this in mind, we will focus on a region of 
parameter space in which we expect to have first order PTs with large supercooling. Supercooling 
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is an example of the hierarchy of mass scales which we argued in Section 4.2 helps to obtain a larger 
CC effect. 



We generalize the Z2XSM potential Eq. (4.8) by relaxing the Z2 symmetry. This allows us 
to write down the three additional operators sh 2 , s 3 , and s, but we eliminate the tadpole by an 



appropriate shift in the field space. We are left with the xSM renormalized potential 



nv>.s\)=${h* 



,.2\2 



1 



1 



(4.34) 



4 2 s 3 2 

The thermal effective potential V e s is calculated in Appendix [F] With this parametrization, 
V e s({h c , s c } ,T = 0) has a minimum at {h c ,s c } = {v,0} where V e g ({v, 0} , T = 0) = and the 
curvatures in the h and s directions are m 2 and m 2 respectively. The Higgs vacuum expectation 
value is fixed by electroweak constraints, but the six r eal numbers {m 2 , m 2 , 64, 63, ax, a 2 } are free 
parameters. 

As in the previous section, we compute the bounce action in order to estimate the PT 
temperature Tp T by solving S^ 3 ' /T ~ 140. This calculation is made more challenging by the 
presence of the additional field direction. To obtain we make the approximation that the PT 
occurs along the trajectory s(h c ) satisfying 

dU({h c ,s c },0) 



ds 







and 



s c (v) = . 



(4.35) 



which reduces the problem back to solving for the bounce in one dimension. In the region of 
parameter space on which we are focused, this approximation gives Tp T to within a few percent 



(see Appendix O for details). Note that the empirical formula Eq. (4.21) cannot be applied here, 



because the effective potential is not well approximated by the form Eq. (4.18). 

We have performed a parameter space scan and searched for a region with large corrections 
to the relic abundance from the CC. In the scan we fix the parameters 63 = —20 GeV, 64 = 
0.2, ai = -25 GeV, and a 2 = 0.2, and we vary m\ £ [(65 GeV) 2 , (170 GeV) 2 ] and m| G 
[(40 GeV) 2 , (140 GeV) 2 ]. In order to connect with the intuition garnered from the single field 



model of Section 4.3 we have mapped the xSM parameter space to a single parameter M . This is 



accomplished by restricting to the trajectory given by Eq. (|4.35|) and defining 

d 



M 



dx 2 



V eS ({h c (x),s c (h c (x))},T = 0) 



(4.36) 



x=0 



where x parametrizes the position along the curve s c (h). The parameter M 2 controls the stability of 
the electroweak preserving vacuum: if M 2 > the symmetric phase remains metastable as T — > 0, 
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Figure 8: The fractional deviation in the relic abundance of a 2 TeV WIMP due to the CC at the xSM 
electroweak PT. The parameter M 2 controls the curvature of the zero temperature effective potential along 
the PT trajectory. For M 2 < large supercooling enhances the CC's effect. For M 2 > the PT does 
not occur, and for M 2 < -2500 GeV 2 the PT occurs before freeze out leading to a suppression of the relic 
abundance shift. 



whereas if M 2 < the symmetric phase becomes perturbatively unstable at some finite temperature 
Tq > 0. In this way, the potential depends on the parameter M 2 in the same way as the parameter 



ao from Section 4.3 We cannot map the xSM parameter space to ao directly because the effective 



potential along the trajectory Eq. (4.35) cannot be expressed in the form of Eq. (4.18). 

In Figured! we have plotted c\e\, given by Eq. (3.46), by projecting onto the M 2 axis and 



choosing mx = 2 TeV. For M 2 < the CC has an O (1) impact on the relic abundance. In this 
region, the supercooling is maximal and Tp T > To = O (few GeV). For smaller values of M 2 , the 
CC effect rapidly decreases and drops below 1% for M 2 < 500 GeV. Therefore, in order for the 
CC to have a significant impact on the relic abundance, the parameters of the scalar sector must 
be tuned into a narrow band where supercooling is large. In Figure |9j we have allowed the WIMP 
mass to decrease to 500 GeV. This change lowers the freeze out temperature, reduces the delay S 
between freeze out and the PT, and therefore increases the CC effect. However, this increase is 
small compared with the amount by which c\€\ varies with M 2 in the M 2 < region. For smaller 
values of M 2 , the PT temperature is higher and for the 500 GeV WIMP, freeze out occurs after 
the PT causing the CC effect to be suppressed. These calculations lead us to the conclusion that 
the optimal region of parameter space is one in which the symmetric phase becomes perturbatively 
unstable at a low temperature and the effective potential is concave at zero temperature. We were 
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Figure 9: This figure shows a subset of Figure [8] as well as the CC effect for a 500 GeV WIMP represented 
by squares. As the WIMP mass is reduced, freeze out occurs at a lower temperature. This increases c\t\ for 
M 2 < where the PT temperature is low, but excludes points M 2 < —100 GeV 2 where freeze out occurs 
after the PT. 



unable to find any points with M 2 > in which the PT completes. 
The following is a benchmark parameter point: 

{di, 6 3 , m h ,m s ,m x } = {-25, -20, 128, 91.1, 2000} GeV, {a 2 , b 4 } = {0.2, 0.2} 



M 



-47.7 GeV 2 {T f , T c , T+ T , T~ T , T } = {107, 70.7, 30.0, 13.7, 12.7} GeV, 



cid = 0.390, 



(69.7 GeV) 4 



(4.37) 



The scalar masses are given by the eigenvalues of Eq. (F.16) which are 



Mh = 141 GeV {0.78,0.22} 
M s = 70.7 GeV {0.22,0.78} 



(4.38) 



with the respective squared eigenvectors indicated to the right. 



5. CONCLUSION 



We have considered a way to probe the hypothesis that the present-day, minute CC energy density 
is the result of a tuning between UV contributions of unspecified origin and IR contributions that 
arise from cosmological PTs. Prior to the electroweak scale PT, the UV contribution would have 
been partially uncancelled leaving an O (M^) energy density. It is possible to probe this energy 
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density with the physics of dark matter freeze out provided that the dark matter mass is greater 
than a few hundred GeV. The dark matter relic abundance is increased due to the effective CC's 
contribution to the Hubble expansion rate during freeze out. 

The notion of how an effective vacuum energy (which is Lorentz invariant in the flat space limit) 
can depend on temperature (which manifestly breaks Lorentz invariance) has been clarified. The 
temperature is an approximation to the mixed vacua, inhomogeneous states whose occupation is 
very probable near the time of the PT. This leads to a spatially averaged equation of state that 
is expressed in terms of an effective vacuum energy density that is somewhere between the false 
and true vacuum energy densities. The true inhomogeneous field configurations may also lead to 
additional dark matter freeze out effects that have not been investigated in this paper. This would 
be an interesting avenue for future investigations. 

To provide a generic prediction associated with the established physics and to provide the com- 
putational details missing in j9], we have analyzed the Standard Model with a 115 GeV Higgs and 
a single WIMP dark matter degree of freedom, assuming that the WIMP interaction effects on the 
dynamics of the PT is negligible. We have found that the CC causes an O (1CT 3 ) fractional increase 
in the relic abundance of a 4 TeV WIMPs. This is typical of non-first order PTs. 

We have also investigated minimal singlet extensions of the SM and searched for parametric 
regimes in which the CC effect on the relic abundance is enhanced. We find that a low temperature, 
first order PT with large supercooling is the optimal scenario for maximizing the CC effect. In this 
limit, the effective CC energy density's contribution to the Hubble expansion rate can be comparable 
to the radiation energy density, and the CC effect can become order one. In the context of a generic 
single field model, we find that reaching this limit requires a tuning of the scalar sector parameters 
and the WIMP mass. Without appropriate tuning, either 1) the PT will not occur at all by thermal 
bubble nucleation, 2) the PT will occur before freeze out (when the dark matter is still in equilibrium 
and the CC effect is suppressed), or 3) the CC effect will not be large. 

As a specific example, we have considered the xSM, an extension of the SM that adds a real 
scalar singlet. In that model, we find that the CC may increase the relic abundance by as much as a 
factor of order few. To maximize the CC effect, the scalar parameters must be tuned into a narrow 
band where fluctuations around the symmetric "vacuum" are slightly tachyonic, which allow for a 
long period of supercooling. The magnitude of the CC effect is relatively insensitive to the WIMP 
mass provided that the latter is sufficiently large such that freeze out begins before the PT occurs. 

The tests of CC fine tuning hypothesis are notoriously rare. In the context of a dark matter 
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probe, it is encouraging that parametric possibilities do exist within simple extensions of the SM. 
It would be interesting to further advance this exploration by computing the dark matter impli- 
cations of modified gravity /self-tuning models and comparing the results with those of this paper. 
Furthermore, it would be interesting to cross correlate other astrophysical tests of those modified 
gravity /self-tuning models with the dark matter predictions made within those models. Note also 
that there are other probes of the cosmological constant during a PT such as gravity wave probes |70] 
that will need more development as the gravity wave spectrum calculational technology improves 

[ZMa]. 



Appendix A: Renormalization Scale 

Any measurable quantity is independent of the renormalization scale. Hence, one should not ex- 
pect that the running of the cosmological constant parameter should affect any physical observable. 
Indeed, the running of the other parameters in the Lagrangian will compensate the running of the 
CC parameter to yield the same Too governing the expansion rate H which can be measured for 
example by a test photon redshift. The renormalization scheme and scale does however determine 
the manner in which radiative corrections play a role. Furthermore, in any practical computations 
involving finite order truncation in ft expansion, there is a renormalization scale dependence to next 
order in the perturbation power unless one is able to explicitly keep exactly the terms of the relevant 
order in ft. 

Given that we are computing homogeneous quantities, one might also naively worry that there 
is a coarse graining requirement down to length scales of H . To see why this is not the case and 
to see what renormalization scales would minimize the radiative correction dependence, consider 
the effective action generating the gravitational equation of motion for the metric g: 

JS eB [g] _ JS EH [g\ \ r)„rhf> iS M[9,<t>] 



where Seh is the Einstein-Hilbert action, the matter field schematically written as 4> satisfies the 
appropriate boundary conditions relevant for the matter distribution, and we assume a renormal- 
ization scale at A. Since we are going to resolve the one-particle thermal states with masses of 
order the freeze out temperature Tf, we should have A > Tf. Semiclassically expanding about the 



classical path </>q on the right hand side of Eq. (A.l), we have 



. r ,4 5tp^(x) S 2 Sj^j [g,<fi] I , 

js es [g] _ e i(S EH {g]+s M [g,M)j\f / D^S^e 2 WW l*=*o (A. 2) 
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where the path integral will have the usual perturbative renormalization. Hence, one can consider 
the physical observables to be defined through 



T^(y) = —^=j^^ ( S M \jg,M -Hn{Af I D A 5(f)exp 



4 u4> 2 {x) 5 2 S M [g,(t>\ , 

2 5^(x) !^ + --- 



(A.3) 

Note that in practice, we are expanding g^ u perturbatively about a homogeneous and isotropic 
FRW background before doing the path integral. Hence, the inhomogeneities can be computed 
using classical perturbation theory and the renormalization scale need not be at A = Hpt even 
though it is at length scales longer than Hp^ for which homogeneity and isotropy are typically a 
good assumptions. 



Appendix B: Derivation of Eq. (3.3) 



Start with the thermally averaged Boltzmann equation for nx(t) 

i^M=-H(nl-nf) (B.l) 

which says that nx tracks the equilibrium number density until freeze out occurs at t = tf . Long 
after freeze out, the equilibrium term can be neglected, and the equation asymptotically approaches 

|(n x a 3 )=-M(n x a 3 ) 2 ^. (B.2) 

One can solve for nx(to) by integrating 

nx(to) = . V 3 ' r ■ (B.3) 

! + »*(*/) (if) %dt(ov)% 

The integral in the denominator accounts for residual annihilations of dark matter particles after 

freeze out. The freeze out time tf is not fundamental but instead an artifact of defining when 

the solution deviates "significantly" from the equilibrium distribution. For temperatures away from 

resonances and thresholds, one can typically parameterize (cry) as 

- T 

(av) = a + b , (B.4) 

m x 



where T is the temperature and mx is the mass of the dark matter. To further reduce Eq. (B.3) 
we apply Eq. (3.25), which implicitly defines tf, and approximate nx(tf) ~ n e ^(tf). Then, the 



denominator of Eq. (B.3) satisfies 



v 1 [ i0 ^,u ffl L^ ~ a+b 2^ 



n x (t f )[^\ / dfcM^«I^ _ 4^ »1 (B.5) 
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for Tf ~ mx/20 the freeze out temperature. Using this approximation we can express the relic 
abundance as 



nx(t ) 



after also applying dt = H d In a 



\na Q /a f ^(a/a/) 



-1 



H 



HI 



(B.6) 



Appendix C: Difference Between Entropy and Energy Degrees of Freedom 



In this appendix, we show that as the universe expands adiabatically during radiation domina- 
tion, the relationship g E (T) = gs(T) hold iff 

dlng E dings 







(C.l) 



cilnT dlnT 

where g E is the effective number of degrees of freedom for the thermal energy density and g$ is the 
effective number of degrees of freedom for the entropy density. We also justify an ansatz that can 
be used to relate g E and gg. 

Assume that the CC energy density is negligible so that p ~ pr, which is the case sufficiently far 



before or after the PT. The entropy and energy densities of a gas are related by Eq. (2.1c), which 
can be written as 

p + P-Ts = (C.2) 
where the pressure P of the gas is given by P(T) = —T(T). The functions g E and gs representing 



the number of relativistic degrees of freedom were defined by Eq. (3.9) and Eq. (3.10) and are 
reproduced here for convenience: 



7T /^ m 4 , 2tH 

p= — q E (T)T and s= 

F 30 y ' 45 



gs(T)T 3 



As the universe expands, energy conservation is enforced by 

d (pa 3 ) + P da 3 = . 



Using Eq. (C.2) and Eq. (C.3) this becomes 



d In g E (T) ^ dlnT ^ = Q 
din a din a g E 



which can be resolved as 



dlnT 
din a 



9S_ 
gE 



1 + 



ldlng E {T) 
4 dlnT 



(C.3) 



(C.4) 



(C.5) 



(C.6) 



39 



Next, impose adiabaticity d(sa 3 )/da = by first using Eq. (C.3) to write 
dln(sa 3 ) 



din a 



dings 
dlnT 



+ 3 



9S_ 
9E 



1 + 



1 dinger) 
4 dlnT 



+ 3. 



and then setting this to zero and solving to find 



9S_ 
9E 



1 + 



1 d\ng E _ 1 d\ng s 
4 dlnT 3 dlnT 



This equation implies g E = gs iff 



dlng E dings 







(C.7) 



(C.8) 



(C.9) 



d In T dlnT 
as claimed. 

To obtain some intuition for this theorem consider the SM electroweak PT. Before the PT, 



the entire spectrum is massless and Eq. (C.9) is satisfied exactly so g E (T) = gs(T) = const for 



T > Tpt- After the PT, we can estimate how much difference between gs and g E is required for 
self-consistency and to justify an intuitive parameterization, by considering a hypothetical situation 
in which one can approximate 



9E/s(T)=g E/s (T i )[T/T i \ 



-12K 



(CIO) 



where K is a constant and T{ is an initial condition temperature. Then, one can solve Eq. (C.8) as 



9S(T) 
9e{T) 



1-3K 
1-AK' 



(C.ll) 



Hence, if < K <C 1, we have a situation in which g E (T) decreases slowly as a function of time 
while satisfying both entropy conservation and gs(T) ~ g E (T). Presumably, K can be viewed as a 
leading term in a Taylor expansion regarding gs/gE- Hence, we will approximate 



g s (T) ^ (1 + K) g E (T) 



(C.12) 



even though we are not necessarily making the assumption of Eq. (C.10) throughout the paper. 



Appendix D: Derivation of Tp T , As, and T(a) 

To find T(a), we start with the temperature before the PT Tp T and impose energy conservation 
to solve for the temperature after the PT Tp T . This allows us to calculate As and 62 in terms of 
A/9 CX . Then, we require the entropy per comoving volume S = sa 3 to be conserved before and after 
the PT to find T(a). 



40 



Assuming that there is a negligible change in a during reheating, we can impose energy conser- 



vation at apT- Using Eqs. (2.6), (3.8), and (3.9), energy conservation can be written as 



2 

^ 9e \Tpr) CTpt) 



30 



(D.l) 



which implicitly defines Tp~ T . This equation can be solved analytically by expanding TX 



PT 



Tp T (1 + At) and linearizing in At along with other small quantities. Using Eq. (3.15) to expand 
gg(T) around g^'(Tf), Eq. (D.l) becomes 

^2 r 



7T 

30 



2 

lim h(apr + e) - h(apr - e) [Tp T f + Ap cx « 4^-g^(T f ) {Tp T f At . (D.2) 

.6 — ^0 J OU 



where we have dropped higher order terms. Using Eq. (3.16), the term in brackets is (7/8)Npx- 
Finally, the equation can be solved for At = Tp T /Tp T — 1 to obtain 



Tp T k. T pT 



1 1 

1 + ^31 + 



^i^{T } )(Tp T Y 



(D.3) 



where €31 is given by Eq. (3.20c). As expected, the energy released Ap cx > controls the reheating 
from Tp T to Tp T . Additionally, the reheating is larger when more particles non-adiabatically 
decouple (larger €31), because the latent heat is distributed over fewer degrees of freedom after the 
PT. 

Next we can calculate the entropy density increase at the PT given by 



V8 = 8®(r+ )-s w (r. 



PT' 



2tt2 
45 



{^(T+ T )(T+ T ) 3 - (Tp T ){Tp T f) 



(D.4) 
(D.5) 



Once again we will linearize in the perturbation by expanding g$ using Eq. (3.15) and writing Tp T 



using Eq. (D.3|). This gives 
As 



2vr 2 ( 




45 




2vr 2 ( 


f 9 { z\T f ) 


45 


{ 9 { s S) (T f ) 



\\mh(apT + e) — h(app — e) + 3At > g$(Tf) (T pT )' 



e 3 i + 3 



1 1 

^31 + 7 



Ap c 



4^(T / )(T- T )4j 



As discussed in Section 



C 
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9 { s\Tf)(Tp T y 



(D.6) 
(D.7) 



we can approximate g^'iTf) ~ g^'(Tf). Then finally As becomes 
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Using Eq. (3.206) and noting T pT apT = TVa/ up to higher order terms, we also obtain 

_1 1 Ap cx 



(D.9) 



Both of these equations illustrate that the entropy increase at the PT is controlled by the amount 
of latent heat released and the number of particles that non-adiabatically decouple. 

Finally we will solve the equation of entropy conservation for T(a). The entropy per comoving 
volume S = s a 3 is conserved excepting the entropy injection at reheating which is assumed to occur 
rapidly at qpt- Entropy conservation may be expressed as 



g s (a) T(a) 3 a 3 = g { s s \T f ) Tfa f + @{a - a PT ) a PT 



/2tt 2 



- 



V 45 



As. 



(D.IO) 



which implicitly defines T(a). To solve for T we use Eq. (3.15) to expand gs(a) then linearize in h 
and As to obtain 



T(a) « Tf 



a f 



1 h(a) _ , .If apr\ 3 

9s (Tf) 3 



As 

" / ) 
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9 ( s S> (T f )Tf 



(D.ll) 



Further expanding h using Eq. (3.16), approximating g$(Tf) ~ gE(Tf), and applying Eq. (3.206) 
we obtain the final expression, 

a f 



T(a) 



L f- 



1 + ^32 f(a) + © (a - a PT ) ~ (e 31 + e 2 ) 



(D.12) 



After the PT, the exotic energy component behaves approximately adiabatically. 



Appendix E: Derivation of PT induced change in the degree of freedom 

We begin with the well-known formula for the energy density of a gas of fermions at temperature 
T with N dynamical degrees of freedom: 

The gas has an effective number of degrees of freedom g& given implicitly by p(T) = |q gE(T) T 4 . 
We can parameterize the decrease in gE due to the decoupling of the fermionic gas by writing 

gE(T)=g E (Tf)- 7 -Nf(a/af) (E.2) 

where 
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The temperature T = T(a) is given by Eq. (D.12) to leading order in the perturbations £j. Since 



/ already multiplies a small term in Eq. (3.16), we need only keep the leading factor in Eq. (D.12) 



which is T = Tf at /a. This lets us write Eq. (E.3 \ as 



f(a) 



8 /30 
7 W 



d 3 p E p 



1 



(a/a f ) 4 



(E.4) 



_e T f + 1 e a f T f + 1. 

Note that /(a) increases from /(ct/) = to /(oo) 1. Due to the exponential temperature 
dependence, the transition to / ~ 1 occurs at T ~ and is smoothly steplike over a time scale 
At w 1///. In this discussion we have assumed .E p = y^p 2 + ?ri^r with mjy constant, that is, we 
neglect any change in the mass of the particle as a function of time. This assumption is valid 
sufficiently far after the PT such that the scalar field expectation value and field-dependent masses 
have approximately stopped varying. 



Appendix F: Thermal Effective Potential Details 

We have calculated the thermal effective potential through one-loop order for each of the models 
in Section [4] Our calculation employs the standard techniques |74H76| . and the case of the Standard 
Model is particularly well documented \49\ 1771 178] . As such, we do not feel the need to reproduce 
the entire calculation here. However, we have chosen to use renormalization schemes which are 
convenient for our calculation, but not standardly employed. Hence, we will use this appendix to 
write down the thermal effective potentials for each of the models in Section [4] and to spell out our 
renormalization conditions. 

In calculating thermal corrections to the scalar effective potential, we do not include contributions 
from the dark matter sector. This is an excellent approximation provided that freeze out occurs 
prior to the phase transition (as we have assumed), such that the dark matter is decoupled from 
the plasma during the phase transition. 



1. Thermal Effective Potential: Standard Model 

1/2 

Let h(x) = y2\H'H\ be the radial component of the SM Higgs field and let h c = (h). In 
calculating the radiative corrections, we need not include the contributions from every field in the 
Standard Model. With regards to the non-thermal corrections, light particles which couple weakly 
to the Higgs can be neglected, and with regards to the thermal corrections, particles which are light 
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and do not decouple during freeze out can be treated as massless. Since we expect that freeze out 
will coincide with the PT at a mass scale of about 100 GeV and that residual annihilations will 
occur down to a mass scale of about 10 GeV, we can neglect particles with a mass below that of the 
bottom quark (i.e., 4.2 GeV). We retain the top quark, bottom quark, physical Higgs, and massive 
gauge bosons 9 which have field dependent masses 

Ml M z lw (K)=pef^ h\ (F.la) 
K(K) = =| (H - J) (F.lb) 

ZV 

where m t = 172.6 GeV, m b = 4.2 GeV, m z = 91.2 GeV, and m w = 80.4 GeV [79]. The non-thermal 
corrections can be expressed as functions of the Coleman- Weinberg potential |74j . Regulating in 
(d = 4 — 2e) spacetime dimensions, the unrenormalized potential is given by 

/If 4 / M 2 3 \ 

where C uv = e~ 1 — je + ln47r and [i is the t'Hooft scale. The thermal corrections can be expressed 
in terms of the bosonic and fermionic thermal functions |75[ [80] 

J B (y) = dx x 2 log (l - e-V^+A =-J2^yK 2 (n^y) (F.3a) 

-* n=l 71 

J F (y) = ^ dx x 2 log (l + e~ V x2+y ) =-J2 -^-V K 2 (nVv) ( F -3b) 

where K 2 (z) is the modified Bessel function of the second kind. Putting the pieces together, the 
Standard Model thermal effective potential (through one-loop order and before renormalization) is 
given by 



V^^T) ~§(kl- V r + \sn + \5m 2 h 2 c + ^ 



12 V cw {M 2 {h c )) - 12 V cw {M 2 (h c )) + 3 V cw (M|(/i c )) + 6 V cw {M^(h c )) + V cw (M 2 (/ 



7T 



2 



- 75.75 T 4 + 



jn4 



90 ' 2tt 2 



-12 J F (M, 2 (/i c ) T- 2 ) - 12 J F (M 2 (h c ) T- 2 ) 



9 



3J B {M 2 z (h c )T- 2 ) +QJ B (M^{h c )T~ 2 ) + J B (M 2 (h c )T- 2 ) \ + O (h 2 ) (FA) 



We work in the Landau gauge (£ = 0) for which the scalar polarization mode and ghost propagators are independent 
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where 8Q, 5m 2 , and 5 A are counterterms. We have also included the term (— |q 75.75 T 4 ), which 
represents the thermal radiative contribution from light quarks, leptons, and massless gauge bosons 
which are relativistic at temperatures T > 10 GeV. The renormalization conditions, 



(h c ,0) 



dh c eff 
—V (SM) (h 0) 



h c =v 



h c =v 







m h 



v£ M \v,o) = o, 



(F.5a) 
(F.5b) 
(F.5c) 



are chosen such that tadpole graphs vanish and \Qg {h c , 0) has a minimum at /i c = v, self-energy 
graphs vanish and the Higgs mass 10 is m^, and the CC is tuned against the vacuum energy density 
to zero. 



2. Thermal Effective Potential: Z 2 xSM 



The Z2XSM potential was specified by Eq. (4.8). Since we focus on the case (s) = 0, we need 



only calculate the effective potential as a function of h c and not s c = (s). That is, the presence of 
the singlet in this model simply add an additional degree of freedom, with field dependent mass 

M 2 s {h c ) = {m 2 s - a 2 v 2 ) + a 2 h 2 c , (F.6) 

to the radiative corrections. We can construct the effective potential from the SM effective potential 



Eq. (FA) as 



K? 2XSM) (hc,T) =v}* M \h c ,T) + V cw (M 2 (h c )) + ^ 2 J B (M 2 (h c )T~ 



rp4 
27T 2 ' 



(F.7) 



An additional UV divergence arises from the term V CW (M 2 ), and is cancelled by solving the renor- 



malization conditions Eq. (F.5) once again for the counterterms. 



3. Thermal Effective Potential: Generic Singlet 



For the theory specified by the action Eq. (4.13), we have the field dependent masses 



M 2 {ip c ) = M 2 - 6£V C + 3A^ C 
M 2 t (Lp c ) = {mi + hnp c f . 



(F.8) 
(F.9) 



10 Since the effective potential is computed from diagrams with zero external momentum, the mass <9^ c T4ff(/i c = 
v, 0) = m\ differs from the Higgs pole mass by logarithmic corrections |81| . which we verify are O (few %). As 
such, we will neglect this distinction and continue to refer to mj as the "Higgs mass." 
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We construct the thermal effective potential as 



V.. 



(GS) 



1 



(<p c , T) = Pex + ^mV' - £<pI + -M + 



cfr 



2vr 2 



N 



j b {mI(^ c ))-aY j Jf{mI i ^ c )) 



(F.10) 



f 1 A\ ^ 1 

+ | ,50 + 5t ip c + ~5M 2 ^ - 88 <p 3 e + —<pi + V cw (M> c )) " 4 X>™ ( M ^)) j 

where 50, <5t, <5M 2 , <5£, and 5 A are counterterms. We do not renormalize using the same renormal- 
ization conditions as we did for the SM. To simply the discussions of Section [4] we have attempted to 
choose the renormalization conditions such that the effective potential preserves certain features of 
the renormalized tree-level potential. For example, the renormalization conditions that we applied 



to the SM, Eq. (-F.5 ), ensured that the effective potential and the tree-level potential agreed to order 
h\ as an expansion around h c = v. In our analysis of Section 4.3, we found it convenient to define 
the parameter «o which controls the shape of the effective potential. This parameter is defined 
using the tree-level potential U((p), but we claim that it also describes the shape of the one-loop 
effective potential provided that the radiative corrections do not significantly distort the shape of 
the potential. For the tuned limit < ao <C 1, this parameter is particularly sensitive to the shape 
of the potential near the origin (p c = since the barrier is very small. The radiative corrections grow 
as (p c — > 0, because the fermions tpi become light, but these logarithmic corrections remain subdom- 



inant. However, with a renormalization scheme of the form of Eq. (-F.5), the counterterms pick up 
a finite piece, which depends on derivatives of logarithms at the renormalization point <p c ~ v, and 
which contributes non-negligibly near (p c ks 0. If we were to use such a renormalization scheme in 
the limit where U((f c ) has a small barrier so < ao <§C 1, then the radiative corrections may lift the 
minimum at (p c ~ and eliminate the barrier. Of course, there is nothing incorrect with using such 
a renormalization scheme except that it is inconvenient since we would not be able to characterize 
the shape of the potential using ao derived from U(tp c )- 

In light of this discussion, we will use a renormalization scheme which preserves the location 
of the minimum at (p c = v and also preserves the shape of the potential near ip c = 0. This is 



accomplished by first writing Eq. (|-F.10|) for T = as 

,(GS) 



V. 



1 



eff (<p c , o) = n(<p c ) + t(vc)<Pc + 2 M (<Pc)<pi - £(vc)<pt + 



Hfc) 4 



(F.ll) 
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where 



N 



16 



8=1 



% c ) = <ft + 



4^2 



A" 



-EM 2 U(vc)-Y, m ^U^ 



1=1 



M 2 (<^ c ) = M 2 + 5M 2 + 
£ (y> c ) = 5 + <5£ + 



ft 

ft r9 



9 



N 



M 2 XU( Vc ) + -£ 2 U(Vc) - 3^Tm 2 ft 2 /^ c ) 



i=l 



N 



4vr 2 

A(^ C ) = A + 5A + A 



i=i 

9 * 



i=l 



and 



In 



In 



M 2 (y c ) 3 
fi 2 2 ' 

Ml^c) 3 



^ 2 

Then the renormalization conditions can be expressed as 



C u 



M 2 (v 





M 2 

£ 

A. 



Near (p c ~ 0, the radiative corrections are at most logarithmic. 



(F.12a) 
(F.12b) 
(F.12c) 
(F.12d) 
(F.12e) 



(F.13) 
(F.14) 



(F.15a) 
(F.15b) 
(F.15c) 
(F.15d) 
(F.15e) 



4. Thermal Effective Potential: xSM 



In Section 4.4 we wrote down the xSM renormalized potential in Eq. (4.34). For general ft, 



and s c , the Higgs and singlet fields mix. In order to calculate the radiative corrections, we must 

h ' 



generalize the field-dependent Higgs mass M?, given by Eq. (F.lb), to the Higgs-singlet mass matrix 
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M? , which has components 



M s ({h c , s c })] u = m 2 /(2v 2 ) (3h 2 c - v 2 ) + s c ( ai + a 2 s c ) 
[M 2 s ({h c ,s c })] l2 = [M 2 s ({h c ,s c })] 2l = h c ( ai + 2a 2Sc ) 
[Ml s ({h c , s c })] 22 =m 2 + a 2 {h 2 c - v 2 ) + 2b 3 s c + 36 4 s 2 c . 



(F.16a) 
(F.16b) 
(F.16c) 



Now we can write down the thermal effective potential in terms of V^g M ^ by subtracting the con- 
tribution from the SM Higgs and adding the contribution from the mixed Higgs and singlet. Doing 
so we obtain 



v}f M \{h c ,s c },T) =v£ M \h c ,T) + b jsi + \m 2 s s 2 c + b jsl + \s c (h 2 c - v 2 ) (a, + a 2 s c 



6b 



5b 3 



1 



V cw {M 2 (h c )) + Tr V cw {M 2 s ({h c , s c })) 



T 1 



+ ^2 t- J 5 (M^h c )T~ 2 ) +Tr J B {M 2 hs ({h Cl s c })T- 2 )} 



(F.17) 



where 5f2, 5bi, and (5aj are counterterms. The trace is interpreted to mean evaluating V cw or Jb with 
f 2 



the eigenvalues of M? . We generalize the SM renormalization conditions Eq. (F.5) to incorporate 



the additional fields, 



dh 



v}f M) ({h c , Sc} ) 



V" / o N "- 

{u,0} V^c/ V^c 



^({^ c ,s c })l { „,o} 



(F.18) 



{n ft , n s } = {1, 0} , {2, 0} , {0, 1} , {0, 2} , {1, 1} , {1, 2} , {0, 3} , {0, 4} , {0, 0} 

where U({h c , s c }) is given by Eq. (4.34). Once again, we require v}f M \{v, 0} , 0) = which tunes 
the CC. 



Appendix G: xSM Bounce Calculation 



As discussed in Section 4.4 the xSM electroweak PT is first order in the parametric regime 
of interest and proceeds through thermal bubble nucleation. In order to determine the bubble 
nucleation temperature Tp T we estimate the action of the three dimensional bounce S^(T) and 
require /T\ T - ~ 140. The bounce field configuration 4>B( r ) is a saddle point solution of the 



PT 



Euclidean equation of motion with an 0(3) symmetry. Let <p = {h, s} be the field space coordinate 
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and let 



■'sym 



v^(T) and 



^brk 



y( b )(T) be the location of the symmetric and broken phases at 



temperature T. In this notation, the field equation and boundary conditions can be written as 

V $ V eG ($,T) = (G.l) 
, lim ${r) = $ sym (G.2) 



d 2 4> 2 dcf) 
dr 2 r dr 



dr 



r=0 



where r is the radial coordinate and V e g is the thermal effective potential. The bounce solution is a 
curve 4>B( r ) which starts nearby to 4>brk at r = and approaches 4> sy m as r — > oo. Once the solution 
4>B{ r ) is obtained, the bounce action is calculated as 



fC 

5 (3) (T) =4vr / 
J o 



r 2 dr 



'>B 



dr 



+ V eS (<p B (r),T) 



(G.3) 



It is difficult to solve Eq. (G.l) by brute force numerics, because the solution is unstable to per- 



turbations about the initial point ^(0), and the over shoot / under shoot method is non-trivial to 
apply in two dimensions. 

Profumo et. al. |59] have outlined a numerical procedure which reduces the calculation to 



iteratively solving the one-dimensional analog of Eq. (G.l). They suggest that one should decompose 



the field equation into a basis with unit vectors parallel and perpendicular to the solution curve (j)(r). 
Suppose that there exists a curve 4>(x) that interpolates between <^>(0) = (f) sym an d <P{L) = <frbrk- Let 

(G.4) 



be the distance along the curve such that L is the total length and 



'II 



dx 



and 



1 



-1 o. 



dx 



(G.5) 



are the unit vectors parallel and perpendicular to the curve at x. In this basis, Eq. (G.l) becomes 

id 2 x 2dx dV($(x)) \ , 



dr 2 ^ r dr 



d 2 4> 



dx 2 



dx\ 
dr ) 



dx 
VrV 







e ± = 0. 



(G.6) 
(G.7) 



The authors of |59] solve these equations numerically using an iterative procedure. 

Since we compute Tp T by calculating the bounce at various temperatures in order to solve 
jrp ^ 140^ ^ ne it era ti ve procedure is too computationally intensive for our purposes. Fortunately, 
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Figure 10: Comparisons of bounce calculations for the xSM benchmark point Eq. (4.37). On the left, the 



bounce action computed at various temperatures between T = 12.7 GeV and T c = 70.7 GeV using the 
method of [59J (squares) and our approximation (circles). On the right, the xSM thermal effective potential 
at Tp T = 13.7 GeV. The solid curve shows the trajectory 4>b(x) obtained using the method of [59 , and 



the dashed curve shows the approximation <pa.pp(x) given by Eq. (G.8). The curves do not coincide at small 



h because the minimum along the h — axis shifts as the temperature is raised. Nevertheless, the action 
along the two paths still agrees remarkably well. 



in the parametric regime of interest the bounce solution 4>B( r ) can be approximated by (fi a pp(x 
{h(x),s(h(x))} where s(h) satisfies 11 

dU({h, s} , 0) 



ds 







and 



s(v) = , 



(G. 



U is the classical potential, and <^ a pp(^) is parametrized by its length x given by Eq. (G.4). Using 
(f> a pp(x), we solve Eq. (G.6) for x(r) and calculate using Eq. (G.3). 



To check our approximation, we also compute the PT temperature using the method of [59] for 



a few parameter sets. In Figure 10 we contrast our approximation with the procedure of [59j for the 



xSM benchmark point Eq. (4.37). We find that our approximation tends to overestimate by a 



few percent generically. However, is a rapidly increasing function of temperature, and even an 
O (5%) deviation in does not causes Tp T to deviate appreciably. 



In the parametric region described in Section 4.4 the solution of dU/ds = is not generally a single- valued 
function of h. However, the boundary condition s(v) — selects out a unique trajectory which tends to stay in the 
"valley" connecting the two minima and passes through the saddle point. 
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